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Abstract. A self-organized model with social percolation process is proposed to describe the propagations 
of information for different trading ways across a social system and the automatic formation of various 
groups within market traders. Based on the market structure of this model, some stylized observations 
of real market can be reproduced, including the slow decay of volatility correlations, and the fat tail 
distribution of price returns which is found to cross over to an exponential-type asymptotic decay in 
different dimensional systems. 
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Recently, microscopic models for financial markets have 
attracted more and more interests @,|,|,|,|,||,|. These 
models are based on the empirical findings of high-frequency 
market data, such as the fat tails in the distribution of 
price changes, the fast decaying of linear correlation, as 
well as the existence of long-term volatility correlations 
@;|Hli§li an d the intrinsic structure of financial markets, 
including the mutual interactions among market partic- 
ipants through herding and imitation behaviors ||[7| or 
switching from one strategy group to another B. 

Although different mechanisms are used in these mi- 
croscopic models to simulate the price formation processes, 
some of the stylized observations of real markets can be re- 
produced. In one set of these models, the Cond-Bouchaud 
herding model j(| and the related percolation models JTo| , 
[Ti"| ] , the power-law asymptotic behavior in the tail of price 
return has been obtained, with an exponent well outside 
the stable Levy regime || as found in real data |l2],[l3|. 
Recently, this Cond-Bouchaud percolation model has been 
modified by introducing a feedback mechanism between 
price return R and trader activity a: a — > a + aR, where a 
is the factor representing the sensitivity to price fluctua- 
tions Jll| . Then the volatility clustering can be produced, 
and more interestingly, the empirically observed asymmet- 
ric pattern of sharp peaks and flat troughs in the market 
prices |i~4| can be incorporated. 

In the percolation-type model |0,|T3,O|, investors are 
randomly distributed on a lattice with certain concentra- 
tion and through the neighboring connections they form 
the percolation clusters, corresponding to herding groups 
or companies with probability a of trading (activity) which 
is set to be the same for all groups. However, the con- 
centration of investors on the lattice is to be pre-setted 
for obtaining resonable results. Here we introduce a self- 
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organized model for information spread based on the social 
percolation process where the investor groups with 

various trading activities and sizes are formed automati- 
cally as a result of the propagations of trading information 
across a social network. 

Consider a d dimensional lattice representing a social 
system with each site located by a trader, that is, the 
whole lattice is occupied. When participating in the mar- 
ket each trader obeys a certain trading way, which can 
be simply represented by the trading activity measuring 
the frequency of putting the buy or sell orders, and she 
may change her behavior due to the influence of outside 
information about trading ways if she believes that the 
information is attractive enough and the new trading way 
can lead to more profits. Suppose that the information 
of various trading ways is not public, and spreads among 
traders according to the communication structure of the 
social system, which for simplification, is considered as the 
transfer along nearest-neighbor links of lattice. Thus, a 
higher dimensional lattice corresponds to the society with 
more connections among agents [i6j . All the traders ac- 
cepting the same information form a cluster or group and 
will exhibit the same behavior in the market as they obey 
the same trading way. 

First, we study the process of information spread for 
different trading ways, characterized simply by activity 
a(r) at each spread step r (and the sensitivity to price 
change a (t) if considering the feedback effect as described 
above In the way of larger a(r) the trading will be 

made more frequently. Each information of trading way 
has an attractiveness q(r), related to trader's personal 
opinion about its profitability. The dynamics of the model 
is: 

— At one step r, a new information of trading way with 
attractiveness q(r) and activity a(r) both determined 
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randomly appears and starts at few informed sites to 
spread by nearest-neighbor connections among the trad' 
with personal preferences pi (both the values of q and 
Pi are set between and 1 and Sp shown below is a 
small fixed parameter). 

— Each step r consists of consecutive substeps of the 
spread between traders: After receiving the informa- 
tion from one of her neighbors at one substep, trader 
i will accept it and obey this new trading way if its 
attractiveness is larger than her personal preference, 
i.e., q(r) > pi. Then she will spread the information 
further to the other neighbors, and increase her pref- 
erence by <5p, i.e., p, — > p, + Sp, as it should be more 
difficult and needs larger attractiveness for someone to 
change her mind at the next step r + 1 to another new 
way. Otherwise, the trader i will keep the old state, 
including the old trading way and preference, and will 
not spread the message to her neighbors until the end 
of this step r. Then, the next substep begins with the 
similar judgement and spread procedure. 

— The spread of one step r will continue until there is 
no source of information for propagation and the pro- 
cedure stops by itself, that is, all the neighbors of all 
the traders accepting the information have already re- 
ceived it, and a trader group or cluster with members 
accepting the common information is formed. Then the 
next step r + 1 starts with the incoming of another 
new information of attractiveness q{r + 1) and activity 
a(r + 1), also selected randomly between and 1, and 
has the spread procedure similar to above. In partic- 
ular, for a trader i belonging to a cluster of previous 
information, when she receives the new information 
with q(r + 1) > pi she will accept it and change her 
trading way to the new one. Consequently, the size of 
that old cluster will decrease by 1 with the expansion 
of a new group. Otherwise, the trader will still remain 
in the old cluster. Thus, if the new information is at- 
tractive enough, the spread of it will form a new group 
which may invade the old ones or even make some of 
them disappearing. Up to now, all traders within the 
same cluster share the same activity a(r) which was 
valid when that cluster was formed. 

Compared with the social percolation model |l5|,[l6|, 
here there is no feedback on the attractiveness q of infor- 
mation, instead, new information with randomly selected 
attractiveness appears at each new step. Moreover, the 
preferences of investors not accepting the new information 
or uninformed are not changed, instead of decreasing, and 
then the trader preferences increase with the acceptance of 
more and more information. After a large number of steps 
for the spread of different trading information, the trader 
preferences of the system are close to the upper limit 1 
and the newly appearing information has little influence 
on the trader structure. Thus, through a self-organized 
process the market traders in this social system automati- 
cally form groups or clusters of various sizes and activities 
due to the acceptance of different trading information, and 
these clusters may correspond to the herding groups of in- 
vestors who imitate each other or individual groups with 



members using the same trading activity a and acting in 
abhe same manner. Here the groups or clusters with larger 
size correspond to the agents in the real market who have 
bigger influence, and different activities of groups repre- 
sent the phenomenon of the real market that the agents 
have different trading timescale and some of them trade 
more frequently than the other. 

Fig. [I] shows the results of our simulations for differ- 
ent hypercubic lattices with length L and dimension d, 
where n s denotes the average number of clusters contain- 
ing s sites, and only two parameters are to be chosen: 
5p = 0.001 and r = 10 5 steps. Initially the trader prefer- 
ences pi are distributed randomly between and 1, and 
each step starts with one randomly selected site on one 
boundary. Thus, the sites close to this boundary represent 
the more informed traders in the market. A Leath-type al- 
gorithm p7| is used for the spread procedures across the 
lattices with helical boundary conditions. From Fig. |l| a 
power law behavior similar to the percolation structure is 
obtained, but only for small and intermediate sizes of clus- 
ters. Compared with percolation theory |l8| ], the effective 
power law exponent in the present model is smaller and 
can be lower than 2 in low dimensions with a finite value 
of ^2 sn s which is very close to 1. This phenomenon for 
the slower decaying of n s may be attributed to the much 
heavier tail in the distribution of cluster sizes in present 
model. In the percolation lattice a cluster with size roughly 
comparable to the total system size (largest or "infinite" 
percolating cluster) appears at or above the percolation 
threshold, while in the model here an earlier percolating 
cluster may be invaded by new-coming information that is 
attractive enough compared with traders' preferences, and 
then its size is diminished with the formation of more and 
smaller clusters. Thus, in the present model the sites be- 
longing to the same cluster need not be nearest-neighbors 
to each other as in percolation models |10[ [, For large s, it 
is not clear whether an asymptotic stretched exponential 
behavior similar to percolation structure above the criti- 
cal threshold exists in present model, and more detailed 
work is needed. 

Next, we study the trading behaviors of these market 
groups or clusters that are determined after enough steps 
of the spreading process as described above. For simplic- 
ity, during the trading and price formation processes sim- 
ulated below, the influence of new-coming information, 
i.e., the change of market structure for trader groups, is 
not considered. The trading of each group j is represented 
by an activity aj obtained from different information in 
the spread process and randomly distributed between two 
limits a max and a m ; n , and then at each time step, group 
j decides to make orders with probability cij (with equal 
probability a.j/2 for buying or selling) or be inactive with 
probability 1 — a 3 . This mechanism does not require that 
each buy order should match a sell order, and the bal- 
ance can be made by assuming the involvement of market 
makers outside the model. Assume that all investors in one 
group contribute the same amount of trading orders, and 
then the trading amount of that group is proportional to 
the group size. Moreover, the price return R at each time 
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step t (used for trading process and different from the 
spread step t), which is defined as the change over time 
interval At of logarithm of pricc^Js supposed to be propor- 
tional to the excess demand HQ, that is, the difference 
between demand and supply obtained from the sum over 
all active clusters, and normalized by the system size L d , 
that is, R — > R/L d , to avoid the influence of lattice size. 

Figs. H and || give the results for the probability dis- 
tribution of price returns on the simple cubic lattice, and 
similar behaviors are found for other dimensions from two 
to seven. In the simulations, the unit of time increment 
At is to be chosen appropriately in order to avoid the 
fluctuations due to random decisions, and here we select 
5 time steps as the unit interval At = 1. Fig. || shows a 
semi-log plot of the probability distribution of normalized 
returns, defined as r — (R — {R))/a with the average of 
return (R) (about 0) over the time series and the volatil- 
ity a = ((R 2 ) - (R) 2 ) 1 / 2 , for a max = 0.02, a min = 0.0001, 
and different time scales At = 1, 2, 10, and 100. With 
the increasing of time interval At, a crossover toward the 
Gaussian distribution is observed from the figure, in agree- 
ment with the finding of empirical financial data fl3|| . 

In our simulations the asymptotic form of the fat tail 
distribution for small time scales is exponential-type, faster 
than the power law behavior with exponent about 4 found 
in recent empirical studies Jl2|,[l3| and theoretical simula- 
tions [Q,|8[|l^]. Similar phenomena can be found in many 
real systems where the probability distributions often cross 
over to an exponential-type decay and exhibit curvature 
in the log-log plots after a limited range of scales for 
power law behavior, and a stretched exponential descrip- 
tion has been shown to account well for many natural 
and economic distributions [^o| . Very recently, an ana- 
lytic form for the whole range of probability distribution 
and the corresponding Langevin equation have been de- 
rived pl| based on a theorem of general stationary ran- 
dom process and the Hong Kong Hang Seng Index data, 
and the asymptotic behavior for large price changes is an 
exponential-type decay: P(r) ~ exp(— /?|r|)/|r|. Fig. || is 
transformed from Fig. ^| by plotting In |r|P(r) as a func- 
tion of the absolute normalized return |r| for the data of 
time interval At = 1. Good fits to straight line in the 
tail region are found with /3 = 0.96 ± 0.02 (positive tail) 
and 0.96 ± 0.03 (negative tail) (fitted over well averaged 
region) . 

However, in above simulations the persistence of volatil- 
ity correlation is absent, similar to the standard Cond- 
Bouchaud model ||. To reproduce the empirical facts of 
volatility clustering and asymmetry of bubbles and crashes 
in stock markets, we follow the method in |ll| to introduce 
a feedback between price fluctuations and group behaviors 
at each time step: aj(t + 1) = a,j(t) + ajR, that is, the in- 
vestors are encouraged by the price increasing with more 
tradings, and are more prudent when price decreases. We 
keep a,j in the interval from a m i n to a max . Furthermore, as 
the trading ways of the groups are different with different 
viewpoints on the price changes, in our simulations the pa- 
rameters etj are set randomly between two limits a max and 
amin- Figs. H to || show the results of simulations in two, 



three, and seven dimensions, where a max = 0.1, a m i„ = 0, 
and the cutoffs of activity a max = 0.2 and a m j n = 0.003. 
To make the initial values of group activities irrelevant, 
the first 10 6 time steps of the simulations are skipped. 

The properties of fat tails in the probability distribu- 
tion of normalized price returns r are shown in Fig. |J, for 
different simulations of time interval At = 1 in square, 
cubic, and d = 7 hypercubic systems, respectively. The 
exponential-type tail behaviors can also be observed in 
the figure, and similarly we transformed Fig. ^ into Fig. 
[| by plotting In rP(r) vs normalized positive returns r, 
with the similar behaviors for negative ones. In the tail 
region good fits to straight line similar to that of the sim- 
ulations without feedback effect are obtained, with pa- 
rameter value (3 as 1.1 for different dimensional systems 
(fi = 1.05 ± 0.01 for size 1001 2 , = 1.108 ± 0.007 for size 
101 3 , and (3 = 1.109 ± 0.007 for size 7 7 , fitted over well 
averaged data). Simulations for other dimensions give the 
similar results. 

The persistence of long-range volatility correlation, which 
is defined as the autocorrelation between the absolute 
value of price return, i.e., [(\R(t) | \R(t+T) | ) - (\R(t) \)(\R(t+ 
T)\)]/[(\R{t)\ 2 } - (\R(t)\) 2 }, is shown in Fig. | for 2, 3, 
and 7 dimensional lattices and time scale At — 1, while 
the linear autocorrelation for the price returns is around 
zero. The inset of Fig. ^| gives the slow decay of volatility 
correlations in log-log plots, representing an asymptotic 
power-law with an exponent about 0.4, similar to the em- 
pirical observation g]. Volatility clustering as well as the 
asymmetry between sharp peaks and shallow valleys of the 
prices found empirically J14| are shown in Fig. [7L which 
gives the time series of one simulation on a 31 simple 
cubic lattice. 

Compared with the Cond-Bouchaud percolation mod- 
els, where the system is fixed at (or slightly above) the per- 
colation threshold or set at varying concentrations with re- 
sults to be integrated over |in,ni|, the model here can au- 
tomatically generate the market structure among traders, 
i.e., the trading groups or clusters, based on the social per- 
colation process which is intrinsically different from the 
Cond-Bouchaud one, and the largest cluster need not be 
omitted for obtaining reasonable results. Moreover, in the 
present model the trading activity for different groups are 
different, reflecting the fact that in real market the agents 
have various timescales or frequencies of trading. 

The qualitative results of the model here are robust 
with respect to parameters. The simulations with differ- 
ent parameters all lead to the same stylized features, in- 
cluding the properties of cluster size distribution after 
large enough steps of information spread, the fat tail of 
price return distribution described by exponential-type 
asymptotic decay, the long-range persistence of volatility 
autocorrelations, and the asymmetry between peaks and 
troughs of the prices. Different dimensional lattices sim- 
ulate the different social societies with various degree of 
personal connections, and the dimensionality of the sys- 
tem does not influence the major results of the model. 
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Fig. 1. Log-log plot of the cluster numbers containing s sites after r = 10 5 steps of the spread model with Sp = 0.001, for 1001 2 
square (averaged over 230 lattices), 101 3 simple cubic (190 lattices), and 7 7 hypercubic (100 lattices) systems. The slope of the 
straight line is —2. 
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Fig. 2. Semi-log plot of the probability distribution of the normalized returns for different time intervals At — 1, 2, 10, and 
100 on 101 3 cubic lattice with a max = 0.02 and a m i n = 0.0001 (averaged over 190 lattices). A crossover toward the Gaussian 
distribution (solid line) is shown with the increasing of time interval. 
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Fig. 3. Replot of Fig. |^ for time interval At = 1 by showing In |r|P(r) as a function of the absolute normalized return \r\, with 
a straight line of slope —0.96. 




Fig. 4. Semi-log plot of the probability distribution of the normalized At = 1 returns r with a max = 0.1, Q^min — 0, ttmax — 0.2, 
and a m i n = 0.003, for 1001 2 square (200 lattices), 101 3 cubic (100 lattices), and 7 7 hypercubic (100 lattices) systems. The 
feedback mechanism between price fluctuation and group activities is considered. 
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Fig. 5. Replot of Fig. H for In rP(r) as a function of the normalized positive return r, with a straight line of slope —1.1. 
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Fig. 6. Autocorrelations for the volatility in different simulations with the same parameters as Fig. [l| The slow decay is 
displayed, and the inset gives the corresponding log-log plots which exhibit an asymptotic power-law behavior with an exponent 
about 0.4 for different systems. 
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Fig. 7. Time series of the logarithm of price in one simulation with the same parameters as Fig. ^, except for a m i n = 0.0003, 
on a 31 cubic lattice. 



